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ABSTRACT 

In addition to its spectrum and temperature anisotropy, the 2.7 K Cosmic 
Microwave Background (CMB) is also expected to exhibit a low level of 
polarization. The spatial power spectrum of the polarization can provide 
details about the formation of structure in the universe as well as its ionization 
history. Here we calculate the magnitude of the CMB polarization in various 
cosmological scenarios, with both an analytic and a numerical method. We 
then outline the fundamental challenges to measuring these signals and focus 
on two of them: achieving adequate sensitivity and removing contamination 
from foreground sources. We describe the design of a ground-based instrument 
(POLAR) that could detect polarization of the CMB at large angular scales in 
the next few years. 

1. Introduction 

The 2.7K Cosmic Microwave Background (CMB) radiation is a vital probe of all 
modern cosmological theories. This radiation provides a "snapshot" of the epoch at which 
radiation and matter decoupled, approximately 300,000 years after the Big Bang, and 
carries the imprint of the ionization history of the universe. This information tightly 
constrains theories of cosmological structure formation. 

The three defining characteristics of this radiation are: its spectrum, spatial anisotropy, 
and polarization. The spectrum and anisotropy of the CMB have both been extensively 
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studied. The COBE FIRAS has determined the blackbody temperature of the CMB to be 
2.728 ± 0.004 K (Fixsen et. al. 1996), and the COBE DMR has detected spatial anisotropy 
of the CMB on 10° scales of AT/T ~ 1.1 x 10" 5 (Bennett et al. 1996). Ground and balloon 
based experiments have also detected anisotropy at smaller scales; see Scott, Silk, & White 
(1995) for a review of these results. However, the polarization of the CMB has received 
comparatively little experimental attention, despite its fundamental nature. The anisotropy 
and polarization depend on the power spectrum of fluctuations as well as the ionization 
history of the universe in different ways. A detection of polarization would complement 
the detections of anisotropy by facilitating the reconstruction of the initial spectrum of 
perturbations as well as the ionization history of the universe. 

The magnitude and spatial distribution of polarization is determined by factors such 
as: the source of the CMB anisotropy, the density parameter f2, the baryon content of the 
universe Vt B , the Hubble constant H, and the ionization history of the universe. CMB 
polarization is uniquely sensitive to the ionization history of the universe, which includes 
the duration of recombination and the epoch of reionization. The detection of, or a further 
constraint on, the polarization of the CMB has the potential to dramatically enhance our 
understanding of the pregalactic evolution of the universe. 

Similar to the CMB anisotropy power spectrum, the polarization power spectrum 
contains information on all angular scales. Large angular scales (larger than ~ 1°) 
correspond to regions on the last scattering surface which were larger than the causal 
horizon. In the absence of reionization, these scales were affected only by the long 
wavelength modes of the primordial power spectrum. This region of the power spectrum 
was measured by the COBE DMR, and establishes the normalization for models of large 
scale structure formation. Similarly, measurements of polarization at large angular scales 
will normalize the entire polarization power spectrum. Because the anticipated signal size 
is small at all angular scales, polarization measurements pose a significant challenge. While 
signals from large angular scales may be weaker than at small scales, the design of a large 
angular scale measurement is comparatively simple and compact, with potentially lower 
susceptibility to sources of systematic error. A detection, or improved upper limit, at large 
angular scales is a natural first step towards probing the polarization power spectrum on all 
angular scales. 

In this paper we review theoretical arguments which suggest that the ratio of 
polarization to anisotropy should be in the range 0.1% to 10%, at large angular scales. 
Existing upper limits on polarization are higher than, or comparable to the measured 
anisotropy level itself (see Table 1). Measurements of anisotropy, by COBE, and other 
experiments on the level of AT/T cmb ~ 1 x 10~ 5 indicate the required level of sensitivity 
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to polarization must be at least AT/T cm b < 1 x 10~ 6 . Thus, to obtain new non-trivial 
information, either a positive detection, or an improved upper limit capable of discriminating 
between different cosmological scenarios, necessitates extremely precise measurements. 

Current detector technology is capable of achieving the required level of sensitivity. 
However, in addition to achieving high sensitivity it is essential to discriminate the 
polarization from systematic effects, such as non-cosmological astrophysical sources of 
polarized radiation. Space-based missions, such as MAP and Planck Surveyor will produce 
full-sky anisotropy maps, and are expected to achieve the required sensitivity level to 
measure polarization as well. The projected sensitivity levels will allow for per-pixel 
detections of anisotropy with signal-to- noise ratios greater than one. The polarization 
maps from these missions, however, are expected to have signal-to-noise ratios less than 
one for each beam-sized pixel, and will be of lower resolution than the anisotropy maps. 
Fortunately, polarization observations are also possible from the ground; as we will 
demonstrate, polarized atmospheric emission is expected to be negligible. 

This paper will concentrate on strategies for a near-term, ground-based polarization 
experiment, Polarization Observations of Large Angular Regions (POLAR), optimized to 
measure CMB polarization at 7° scales, for ~ 36 pixels. The design incorporates many 
techniques developed for previous anisotropy and polarization experiments, from the 
ground, balloons, and space. The primary goal of the paper is to describe the feasibility of 
measuring large angular scale polarization, and highlight the conclusions which could be 
drawn from such a measurement. In section |] we review the theory of CMB polarization 
which motivates the experimental design. We describe the main experimental challenges in 
section || and focus on two which affect the global design of the instrument: discrimination 
of CMB polarization from polarized foreground sources (in sections f| and |5|), and an 
observing strategy designed to minimize the time required to detect a cosmological signal 
(section H). Finally, we estimate the polarization signal we expect in several different 
cosmological scenarios, and speculate on the conclusions which could be drawn from such a 
detection. 



2. CMB Polarization: Theory 

Anisotropy of the CMB is generated by metric perturbations of the universe. There 
are two primary types of perturbation that generate anisotropy of the CMB: scalar 
contributions, generated by matter and radiation density inhomogenities, and tensor 
contributions, associated with gravitational waves. Both types of perturbation give rise to 
temperature fluctuations in the CMB via the Sachs- Wolfe effect (Sachs & Wolfe 1967). 
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Thomson scattering of anisotropic radiation by free electrons inevitably generates 
polarization (Chandrasekhar 1960). Scattering by a single electron produces polarized 
radiation with an intensity approximately 10% of the anisotropy quadrupole amplitude 
when averaged over all directions of photon incidence and scattering. In the case of CMB 
polarization, the exact polarization level, as well as the angular scale of the distribution 
of polarization on the sky depend on the optical depth along the observer's line of sight, 
and on the particular sources of metric perturbation (Rees 1968; Basko & Polnarev 1980; 
Negroponte & Silk 1980; Tolman 1985). For a recent review see: Hu (1996), Kosowsky 
(1996). 

According to the standard model of the evolution of the pre-galactic medium 
after recombination, the previously ionized plasma formed neutral hydrogen which was 
transparent to the CMB. However, the universe may have undergone a secondary ionization 
of the recombined hydrogen. Gunn & Peterson (1965) formulate a measurement of the 
ionization fraction of the intergalactic medium using the lack of a Lyman-a trough in the 
observed spectra of distant quasars. Recent results show that the majority of intergalactic 
hydrogen is highly ionized to a redshift of at least z ~ 5, indicating that the universe 
must have reionized at an earlier epoch (Peebles 1993). Several models predict reionization 
occurred in the redshift range approximately 30 < z x \ < 70 (Ozernoi & Chernomordik 1976; 
Gooding et al. 1991; Durrer 1993; Tegmark & Silk 1993; Nasel'skii & Polnarev 1987). 

In contrast to the standard model of recombination, non-standard models invoke 
additional non- equilibrium sources of ionization. These models predict a prolonged, or 
even non-existent, recombination and/or subsequent ionization of the recombined plasma. 
Since polarization is generated by scattering of photons on free electrons, its magnitude 
and spatial distribution could be used to discriminate between non-standard models and 
the standard model (Bond & Efstathiou 1984,1987; Basko & Polnarev 1980; Nasel'skii & 
Polnarev 1987; Ng & Ng 1996; Zaldarriaga & Harrari 1995; Crittenden, Davis, & Steinhardt 
1993; Frewin, Polnarev, & Coles 1994). An early reionization effectively introduces an 
additional 'last' scattering surface. This has two effects, both of which, in principle, can 
enhance the magnitude of the polarization on large angular scales. Primarily, the additional 
scattering of photons during reionization can create new, or amplify existing polarized 
radiation via the Thomson mechanism discussed above. Additionally, the second 'last' 
scattering surface occurs at a much lower redshift, implying that the causal horizon on this 
rescattering surface is larger, and thus, will subtend a larger angle on the sky today. 

In general, large-scale polarization is enhanced in models which predict early 
reionization. As we will demonstrate, for reasonable non-standard models, the amplitude 
of polarization on 10° angular scales is on the level of 10% of the anisotropy, while for the 
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standard model of recombination the corresponding polarization level does not exceed 1%. 
It is worth mentioning that all of these models predict approximately the same level of 
anisotropy at 10° scales, and hence all of them are compatible with the results of the COBE 
DMR experiment. 

In the remainder of this section we will illustrate the important theoretical features of 
the polarization of the CMB. We will first describe an analytic treatment which predicts the 
level of polarization for both standard and non-standard reionization h istories. In section 
2.3 we will describe a numerical simulation of the effect of a non-standard reionization 
history on the polarization of the CMB. We will find that the more qualitative analytic 
results agree quite well with the quantitative results of our numerical simulations. 



2.1. Polarization produced by cosmological perturbations 

Here we develop the mathematical formalism which will allow us to describe the 
polarization of the CMB in a consistent fashion. With these tools, we will subsequently 
determine the polarization signal we expect to observe, using two different techniques. 
The first method is an analytic approach which will provide a physical framework for 
understanding the polarization of the CMB. The second approach is more quantitative, and 
will allow us to obtain numerical estimates of the polarization signal. In order to describe 
the polarization of the CMB, we will first introduce a parameterization which describes 
the polarization state of arbitrary radiation fields. We then apply this formalism to the 
polarization state of the cosmological signal which we are seeking to detect. 

Consider a polarized electromagnetic wave with angular frequency, u: 

— * 

E = E y0 sm(ujt — S y )y + E x0 sin(a;t — 5 x )x. The polarization state of electromagnetic 
radiation can be characterized by the Stokes parameters: /, Q, U, and V. I — I y + I x , with 
I y = (Ey Q ) and I x = (E^q). I is the total intensity of the radiation, and is always positive. 
The parameters Q = I y — I x and U = 2E y0 E x0 cos(5 y — 5 X ) quantify the linear polarization 
of the wave, and V quantifies the degree of circular polarization (when V = 0, the radiation 

is linearly polarized or unpolarized). The level of polarization is defined as II = tE + — ) 
and the polarized intensity is I po i = II x 7. 

An alternate representation for the Stokes parameters will be of use in the following 
sections. We introduce a symbolic vector for the distribution function of occupation 
numbers of polarized radiation: n = where I is the symbolic vector introduced in 
Chandrasekhar (1960) and related to the Stokes parameters in the following way: 
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Since Thomson scattering cannot produce circular polarization, V — 0, we will consider 

( h\ 

. An unpolarized distribution in zero-th order approximation is 



the 3-vector: I 



given by n = n 
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As shown in (Basko & Polnarev 1980) and further discussed in (Polnarev 1985; 
Zaldarriaga & Harrari 1995), polarized radiation in the presence of cosmological 
perturbations can be represented as 











h = n Q 
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+ ni 



where fix = riA + fin is the correction to the uniform, isotropic, and unpolarized 
radiation described by n . The Planck spectrum, n Q , depends only on frequency, and 
iia + An ar e the anisotropic and polarized components, respectively, which are functions of 
the conformal time, r], comoving spatial coordinates, x a , photon frequency, u, and photon 
propagation direction specified by the unit vector e(9, <fi) with polar angle, 9, and azimuthal 
angle, 0, given in an arbitrarily oriented spherical coordinate system. 

The Equation of Radiative Transfer in terms of 11(77, x °i v i 0); where \i = cos 9, is: 



<9n a <9n 



dndu - 
dx a dv drj ^ U 



and 



1 r 2n 



An 7-i J 



P(/x, 0, n', <p')n(r), x a , v, n', <j>')dn' dtf , 



(2) 
(3) 



where q = axN e a, and the Einstein summation convention is implied. In these 
expressions, a is the cosmological scale factor, P is the scattering matrix described by 
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Chandrasekhar (1960), ot is the Thomson cross section, and N e is the comoving number 
density of free electrons. In general, the effects of a particular choice of metric perturbation 
are manifest in the first term on the right hand side of (|2[): 

dv_ _ I dhap 
dt] 2 dr) 

(Sachs & Wolfe 1967). After retaining terms up to first order in metric perturbations, h a p, 
and since |^ is of the first order, we can replace |^ by in the source term (v a is the 
unperturbed frequency). This implies that the factor 



7 



u a dn Q d In n c 



n n dv n d In v n 



gives a universal frequency dependence for anisotropy and polarization effects, independent 
of the type of metric perturbations (Basko & Polnarev 1980). 

The angular dependence of P is such that 



(where is the symbolic zero- vector) , so we conclude in the zero-th order approximation, 
J = 6. For the first order approximation, in the following, we will understand by J actually 
Ji, in which n is replaced by h^. 

After linearization and spatial Fourier transformation, the equation of transfer takes 
the following form (with v replaced by v): 



^+^n 1 £ = 7#S-g(n 1 £- J £ ). (5) 

Here, = —^h^^e 13 , and "' " = j-. We have specified spherical coordinates in such 
a way that fi = cos#, where 9 is the angle between a vector along the line of sight, e, and 
the wave vector k and <j) is the azimuthal angle of the vector e, in the plane perpendicular 
to the vector k. 

For a given k, h af3 ^ can be represented as a superposition of scalar waves (below we will 
use subscript "S") and tensor gravitational waves (subscript "T" ). Taking into account 
the tensorial structure of the waves, and restricting our consideration to perturbations with 
wavelengths longer than the cosmological horizon at the moment of equipartion (i.e. at the 
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moment when the energy density of matter equals that of radiation, see for example Harrari 
& Zaldarriaga (1993)), we can write 



H: 



1 



r/k n Ks(k) 



'1 — \x ) cos! 



d 



Id/ sin krj 
r\ drj I ?7 



«t(A;) 



(6) 



Here, the y |^5,t(^)| 2 are the amplitudes of the corresponding metric perturbations at the 
moment when their wavelengths are equal to the cosmological horizon. 

For krj <C 1 we have: 



1 3 
H% ~ TE^V^O) - -(1 - /x 2 ) cos20k t (A;) 



(7) 



while for krj 3> 1, 



1 3 
#n = — rik 2 n 2 Ks(k) + t— ^(1 - /i 2 ) cos 20 cos kr)K T (k) (8) 

15 CT^ 

For a plane wave perturbation with wavevector k, the anisotropy and polarization can 
be described as (Basko & Polnarev 1980): 



n A = «s (a* 2 



1 

V / 



+f(l-/^ 2 ) 



(1\ 

1 

V / 



cos 20 



(9) 



nn = /3s(l - /i 2 



( 1 \ 



V / 



+ (3i 



( (l + /i 2 )cos20 \ 

-(1 + /i 2 ) cos 20 
i 4/isin20 J 



(10) 



Substituting (Q) and ( |TU| ) into the integro-differential Equation of Radiative Transfer, 
we obtain the following system of coupled ordinary differential equations for ols,t and 



S,T- 



3 1 



(11) 



£,S,T + qt,s,T — Fs,t, 



(12) 



- 9- 



where £$,t — a s,r + Ps,t> an d F ST is the appropriate source function for scalar or 
tensor perturbations. This system of coupled equations illustrates the intimate relation 
between anisotropy and the generation of polarization. Integrating this system of equations 
we obtain the following general solution for Ps,t'- 

Ps,t = I f F SjT \e~ T - e-& T l drf (13) 
7 Jo J 

where t(t),t)') = p, q(x a )dx a is the optical depth with respect to Thomson scattering. 

For wavelengths large in comparison with the cosmological horizon at the moment of 
decoupling, r] D , (kr] D <C 1), the source function, at this moment can be approximated by: 

F »-&"{J&> (14) 

It can be shown that the source functions are rather insensitive to the exact functional 
form of the variation of the optical depth with respect to time. (Basko & Polnarev 1980; 
Nasel'skii & Polnarev 1987) These functions are primarily characterized by the epoch and 
duration of decoupling. Following (Zaldarriaga & Harrari 1995), we adopt the following 
approximation for the time variation of the optical depth: 

dr] 

dr = — : T 



Ar] 



D 



(see also (Basko & Polnarev 1980; Nasel'skii & Polnarev 1987) for a more detailed 
discussion.). Here Ar^ is the characteristic time scale of the duration of decoupling. 
Approximating the source functions under the integral ([13|) , by their values at the moment 
of decoupling t]d, which gives the main contribution to polarization, we have 



1 
7 

The integral in (|15D can be evaluated in the following way: 



dr 

Ps,t ^ ^(F S:T )\ D A VD I e- T - e-^ T — (15) 
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Hence, 



n 



1 , 10 A ,2 

105 3 ' ' 





f 1 ^ 


3 I 


J^)(l-^ 2 ) 


-1 

V J 


V 



— (1 + fi 2 ) COS I 

4/i sin 20 



(17) 



Comparing equation ( |TTD with equation (|H]), we find that the polarization generated 
by a single perturbation mode with wavevector k is given by: 



2 10 f or 

% = -Yog In y (VDk)(A D k)j U s (k)(l - /i 2 ) - -kt(A;) [(1 + /i 2 ) cos 20 + 2/isin20 

' (18) 

Now we can calculate the root mean square (rms) polarization measured by an antenna 
with an effective averaging angle A - The main contribution to the rms polarization, n(B A ) 
, is contributed by modes with k < /c max (0A) ~ = 



n(e A ) = J(n 2 ) fe> ^ (19) 



^ In YVDATjDiy/QsBg + Q T B T , (20) 



where: 



B S 



[ (l + /i 2 ) 2 / * cos 2 2<pd<p + [ Vci/x / ^ sin 2 20 
J-l Jo J-l Jo 



9 

8^ 



36 



(21) 
(22) 



Here Q s ,t = fo ' 



k*\K s , T (k)\ 2 f, With \K S>T {k)\< 



K 



S.T 



k ns ' T , and 



arc 



the amplitudes of perturbations with wavelengths equal to the cosmological horizon at the 
present moment (n = corresponds to a scale invariant Harrison- Zeldovich spectrum). 
These amplitudes are normalized to the COBE DMR anisotropy quadrupole detection 
which is approximately equal to 2 x 10~ 5 . Assuming that ns = = n, we obtain 



\ 8 10 . 1 

n(6 A ) = — r— In— 7] D Ari D -=== 
105V15 3 y/4: + n 



360 



O I Tl 

°\ 2 +2 



7 V< + T^ 



(23) 
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Taking into account the relationship between redshift and conformal time, z ~ ^ , we 
have that ^ lAm hence 

ZD Vd ' 

2 Ar? D 1Az d 1Az d 1 /zsd\ 
r] D 2 z D 2 2 D z 5D V Zd ) 

where zsd is the redshift of decoupling predicted by the standard model of recombination. 
Finally: 

n(e.) = 4 x io- 7 ^ (pi \k, 9 (24) 

where 



2 x 10- 2 lnf ^ 7° \ 2 V< + < 10 3 ^ , n\- 1/2 ^360° 



(25) 



105^ l360°J 2xl0- 5 z SD { 1+ 4 J I 7° J ^ 1 + # 2 
and 

9 = 

is the ratio of the tensor perturbation amplitudes to the scalar amplitudes. 

The factor incorporates the perturbation amplitudes, normalized to the anisotropy 
quadrupole measured by the COBE DMR. It contains all information about the type of 
metric perturbation, allowing us to isolate factors which depend upon the nature of the 
perturbations, and those which do not. For n = 0, and g = (i.e., no tensor perturbations), 
K Dj9 ~ 1. When g = oo (i.e., no scalar perturbations), K 0i9 ~ 1.84. Finally, when g = 1 (i.e., 
equal tensor and scalar contributions), we find K Di5 ~ 1.47. From this we observe that 
is rather insensitive to the ratio of tensor to scalar amplitudes, g. 

We now emphasize the angular regions to which the preceding discussion is relevant. 
Equations (^) - (|25|) ( which are based on asymptotic formula (|7]), and the approximations 
used in (|i~5D), are valid for modes which satisfy: kArju < 1. In terms of angle on the sky, 

360° A VD 

We can apply equations (p3|) - ( ^5f) to an observation which has an angular resolution 
0^4, as long as: 

e. > = 360 = (— ) {— ) -6— (— j \- w .) . 
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As an example, the standard model of recombination predicts ~ 0.1, which implies 
that ®A min — 0.6°. For pure scalar perturbations (n = 0), the expected level of polarization 
at this angular scale is: 11(0. 6°) ~ 6 x 10~ 6 . For an observation with 0^ ~ 6°, the 
polarization is n(6°) ~ 5 x 10~ 8 . The observed polarization is suppressed by a factor of 
~ 100 with this lower resolution beam. 



Consider another example, a non-standard model for which =jj^- ~ 1, and zo ~ z$d-, the 
angular scale is: BA mi „ — 6°. The polarization predicted in this scenario is: 11(6°) ~ 5x 10~ 7 . 
Finally, for 9^ < 0yi min , the polarization is suppressed, and its dependence on 9^ is 
determined by the details of the ionization history (Zaldarriaga 1997; Polnarev 1985; 
Nasel'skii & Polnarev 1987; Bond '& Efstathiou 1984,1987). 

To summarize, for a given O^, the polarization level is proportional to (see 



equation fl24|)) and is smallest for the standard model of recombination. Alternatively, this 
analytic approximation applies to smaller angles in the standard model, as opposed to 
the larger angles predicted by non-standard models (see Nasel'skii & Polnarev (1987) for 
a more detailed discussion. Figure [I] schematically illustrates the angular dependence of 
polarization in standard and non-standard models. 



2.2. Polarization Power Spectrum 

The analytic treatment above describes the essential physics responsible for the 
generation of CMB polarization. We have discussed the aspects of non-standard 
recombination which are relevant to the large scale polarization of the CMB. In order to 
estimate the observable polarization signature, we now detail a more quantitative approach 
based on the polarization power spectrum. This approach also allows us to discuss the 
effect of an early reionization on the observed polarization. 

For quantitative estimates, the polarization and anisotropy source terms which appear 
in the equation of transfer can be decomposed into Legendre series. The individual modes 
are then evolved to the present where the spatial structure of the CMB can be computed 
(see Bond & Efstathiou 1984,1987; Ng & Ng 1995, Zaldarriaga & Harrari 1995; Frewin, 
Polnarev, & Coles 1994, for example). Because the CMB is an imprint of the epoch of 
linear evolution of perturbations, the individual modes evolve more or less independently. 
This treatment lends itself particularly well to numerical analysis (Seljak & Zaldarriaga 
1996a). The relevant results of such analysis to the present discussion are the anisotropy 
and polarization power spectra. Here we connect the results of these numerical procedures 
with the analytic treatment presented in the previous subsection. 
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The temperature of the CMB, being a scalar valued function, can be expanded in a 
spherical harmonic series on the sky, at a particular point on the sky, x: 

T((x) = 5^ ar ttm Y im {±) (26) 

£,m 

where the Y^ m (x) are the spherical harmonics at x. The temperature two-point correlation 
function is given by: 

C T,£ = 2£ + \ *52( a T,em a T,£m) (27) 

The variance of the ar,i m is given by the Ct ; i, since Var[a^ m ] = (\aT,£ m \ 2 ) — {\^T,im\) 2 = 
{\a>T,e m \ 2 ) = Cr,e if the ax/m are gaussian distributed with zero mean, and (. . .) denotes a 
whole-sky average followed by an average over all observational positions. 

The polarization of the CMB is a tensor-valued function, with a symmetry group 
different from that of the anisotropy. As shown in (Zaldarriaga & Seljak 1997), complex, 
linear combinations of the Stokes parameters transform under rotation about the line of 
sight by an angle ip as: 

(Q ± zE/)'(x) = exp T2ilp (Q ± iU)(pt) (28) 
The analogous expressions to equation |26| for the Stokes parameters are: 



(Q + iU)(±) = E^m2^m(x) (29) 
i,m 

(Q-iV)(±) = ^^^^(X) (30) 



Lm 



where the spin-weighted spherical harmonics, -^Xtmi are a complete and orthonormal 
set of basis functions on the sphere (for an equivalent technique, which uses second rank 
tensors on the sphere to describe the Stokes parameters, see (Kamionkowski et al. 1996)). 
Taking complex, linear combinations of the expansion coefficients, a± 2em , one defines: 

ClE,em = -(lym + a -2,feJ/ 2 (31) 
O-Bjtm = i( a f,em — a -2,£m)/2- (32) 

From these we construct two independent correlation functions which characterize the 
polarization: 

^E,e - 2l+\ ^^ a *E/m a E,lm) (33) 



2£ + 

Cb,£ = 2£+l ^2( a *B,em a B,£m) (34) 
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The Cgg and Cg e have different physical origins and have interesting properties under 
symmetry transformations, such as parity inversion. We refer to (Seljak & Zaldarriaga 
1996b) for further discussion of these fascinating spectra. 

We now define Cf = e + C# t and form: 

1 oo 

(Q(x!)g(x 2 ) + C/(x 1 )[/(x 2 )} = — £(2^ + l)C?P e (cos6) (35) 

where xi and x 2 are vectors toward two different locations on the sky separated by an angle 
6. When 9 = we have for the polarization autocorrelation: 



/ P 2 oi = (Q 2 + f/ 2 ) = ^E(2^ + i)Q n 



We can now connect the results of this power spectrum calculation with the analytic 
treatment presented above. For a measurement of the polarization at an angular scale A , 
we have for the observable level of polarization: 



n(e A ) = ^ = \ 



1 oo 

- £(2* + 1)< A x <T, (36) 

1=2 



I I\ 4tt 

where I ~ 2.7K, and the window function, W® A } is a factor which weights the contribution 
of the £ th moment to the power spectrum. It quantifies the angular sensitivity of a given 
experiment (see section |6] for further detail). For a given theoretical model (i.e. Cg), 
equation |36] can be compared to the analytic expression given in equation |4j. 

The polarization power spectrum, Cf, is highly dependent on the cosmo logical details 
of the model used to generate the anisotropy spectrum, Cf. For this reason, the ratio 
of polarization to anisotropy is often calculated in order to predict observable levels of 
polarization for particular observations (Ng & Ng 1995; Crittenden, Davis, & Steinhardt 
1993). In practice, the polarization spectra are normalized to the appropriate anisotropy 
spectra, which are in turn normalized to the low I values of Gf measured by the COBE 
DMR. The power spectra are generated numerically by CMBFAST (Seljak & Zaldarriaga 
1996a), which permits simultaneous calculation of anisotropy and polarization, as well as 
consistent normalization. 



2.3. The Effect of Reionization on the Polarization Power Spectrum 



As mentioned earlier, non-standard models of the ionization history are characterized 
by non- instantaneous decoupling and/or non-zero optical depth along CMB photon 
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trajectories. We have discussed the effect a of non- instantaneous recombination using the 
analytic method treated above. We now wish to examine the effect of reionization on the 
details of the polarization of the CMB. This investigation lends itself particularly well to 
the numerical evaluation of the polarization power spectrum, calculated using numerical 
routines such as CMBFAST. 

In general, models of reionization often rely on structures such as an early generation 
of stars (Population III), or energetic proto-galaxies to provide either ionizing radiation or 
collisional heating mechanisms. Thus, for every model of reionization there corresponds a 
structure formation scenario, as well as a commensurate set of cosmological parameters to 
be confronted with observational evidence. We will not speculate here on the plausibility 
of specific models of reionization. Discussion of mechanisms for early reionization can be 
found in (Ozernoi & Chenomordik 1976; Gooding et al. 1991; Durrer 1993; Tegmark & 
Silk 1993; Nasel'skii & Polnarev 1987). As noted above, the Gunn-Peterson test provides 
definitive evidence for an ionized intergalactic medium out to a redshift of at least z — 5. 
In fact, the upper limit on the redshift of reionization is set only by the paucity of observed 
quasars beyond z = 5 and, in principle, could be much higher than this. The COBE FIRAS 
limit on the Compton-?/ parameter y — J drkb(T e — T cmh )/m e c 2 < 2.5 x 1CT 5 (Fixsen et 
al. 1996), severely restricts the energy input allowed in models of reionization, but does 
not tightly constrain the epoch of reionization or the ionized fraction of the intergalactic 
medium. The limit is compatible with many early reionization scenarios. 

The effect of reionization can be parameterized in two equivalent forms. One method 
is specified by the optical depth, r ri , for photons due to Thomson scattering along a line of 
sight to the last scattering surface. The second method specifies the redshift of reionization, 
z T i, and the fractional ionization x (electron-to-proton ratio). The two parameterizations 
are related as follows (Peebles 1993): 

— (t)^(?)"(4)—' 3/2 <"> 

where h is the Hubble parameter, Q is the total density parameter of the universe, and 
CIb is the density parameter of baryonic matter. Equation ( |3"7j ) shows the effect of curvature 
of the universe on the optical depth. For reionization occurring at the same redshift and 
ionization fraction, in an open universe (Q < 1), the optical depth will be greater than in 
a flat or closed universe. We also note that the physical size of regions which are in causal 
contact (Hubble radius) at the epoch of reionization, t ri , is of order ~ ct Ti . We expect 
that regions of this size will produce coherent polarization of the CMB, and affect the 
observed polarization power spectrum at angular scales which correspond to the angular 
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scale subtended by the horizon size at the epoch of reionization. This argument is similar to 
those which predict a coherence scale in the CMB anisotropy power spectrum. For example, 
the acoustic peaks in the CMB anisotropy power spectrum arise from causal mechanisms 
(i.e., sound waves propagating in the photon-baryon fluid) acting on scales of order the 
horizon size at the epoch of decoupling. A similar effect occurs for the CMB polarization 
power spectrum, though in this case it is the horizon size of the re-scattering surface, not 
the 'primary' scattering surface, which is imprinted in the observed power spectrum. 

Following Peebles (1993), we expect that the observed CMB polarization angular 
correlation scale will be: B r i ~ O.l^sO/i) 1 / 3 rad. For O = 0.1, Ob = 0.1, h = 1 we find 
ri ~ 1°, and for Q = 1,Qb = 0.05, h = 0.65 we find 6 ri ~ 2°. This new angular scale, 
absent in non-reionized models, is manifested in the spatial polarization correlation function 
and creates a peak in the reionized polarization power spectra at i < 20. 

Using CMBFAST, we have generated polarization spectra created by 
scalar perturbations in a CDM dominated, completely reionized, universe with 
x — 1,Q — 1,0_b = 0.05, h = 0.65. By varying the redshift of reionization in the range 
< z ri < 105, we compute multiple polarization power spectra, which are displayed in figure 
^|. The power spectra illustrate the main features expected from the theoretical principles 
detailed above. Large angular scales correspond to modes with wavelengths greater than 
the width of the last scattering surface. Prior to recombination photons and baryons were 
tightly coupled and the relatively short timescale for acoustic oscillations prevented the 
formation of long-wavelength perturbations. These effects are particularly evident in models 
without reionization. 

In models with early reionization, polarization at large angular scales is enhanced due 
to multiple photon scatterings following reionization. At smaller angular scales (£ ~ 100), 
in models with and without reionization, the polarization power spectra exhibit oscillatory 
behavior, caused by the same type of acoustic oscillations which generate the Doppler 
peaks in the anisotropy power spectra (Frewin, Polnarev, & Coles 1994; Zaldarriaga 1997). 
Though not relevant for the large angular scale considerations discussed here, for t 3> 100 
the polarization is highly suppressed due to Silk Damping (Hu 1996). 

The power spectra are, effectively, predictions of the polarization which should be 
observable given a particular observing strategy. We will show in section ^| that the 
rms polarization expected from the spectra shown in Figure 2 with = 7°, is in the 
range 0.05/iK< J po i < 1.0/iK, where the lower limit is standard recombination with no 
reionization, and the upper limit is for total reionization starting at z = 105. These limits 
agree well with the analytic estimates for non-standard ionization histories discussed in 
subsection 2.1. For a 6° experiment and a non-standard ionization history, figure 1 predicts 
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a polarization level of 5 x 10 7 ~ 1/iK which agrees well with our numerical simulations of 
early reionization (e.g., for z reionization = 105). 

3. Experimental Overview 

Measurement of polarization of the CMB poses a wide variety of experimental 
challenges, many of which are familiar from the experiments now measuring spatial 
anisotropy in the CMB. We describe below the design of POLAR to illustrate the 
experimental issues that must be addressed in any CMB polarization observation. 

POLAR will measure polarization on 7° scales with two separate radiometers, one in 
the K a frequency band, and one in the Q band, covering the spectrum between 26 and 46 
GHz. These radiometers operate simultaneously, and their frequency bands are multiplexed 
into several sub-bands to allow for discrimination against foreground sources. Each 
radiometer executes a drift scan of the zenith with a separate FWHM=7° beam produced 
by a corrugated feed horn antenna. POLAR will observe ~ 36 different pixels for many 
months to reach the level of a few /iK per pixel. The design builds on techniques developed 
in previous searches for CMB polarization (Nanos 1979; Lubin 1980; Lubin & Smoot 1981; 
Lubin, Melese, & Smoot 1983; Wollack et al. 1993; Netterfield et al. 1995) and is driven by 
the key issues identified in this paper: the size and angular scale of the anticipated CMB 
signals, spectral removal of foreground sources, optimization of the observing scheme, and 
anticipated systematic effects. 

3.1. The Polarimeter 

Radiation from the sky couples into a corrugated circular horn antenna (See Fig. 
|3|). This antenna has extremely low sidelobes, near —80 dB at 90° off axis, in both 
polarizations, across a full waveguide band. The main lobe of the antenna is gaussian with 
a FWHM of 7° which is near the minimum that can easily be obtained without additional 
optical components such as lenses or primary reflectors. The antenna output couples to an 
ortho-mode transducer (OMT), a waveguide device that decomposes the incoming wave 
into two orthogonal linear polarization components. The OMT defines the x-y coordinate 
system of the antenna. 

The Q and U Stokes parameters are defined in terms of a coordinate system fixed to 
the sky. There are several approaches to measuring Q and U for a particular pixel on the 
sky. Lubin & Smoot (1981) employ a Dicke switch which alternately couples each of the 
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polarization components from the OMT to a low-noise amplifier and square-law detector 
(Lubin f980; Lubin & Smoot 1981). Phase- sensitive detection at the modulation frequency 
of the switch yields the difference between these two components, the Q Stokes parameter, 
and helps overcome 1/ / noise from the amplifier. One can show that after a 45° rotation 
about the antenna symmetry axis the instrument measures the U parameter. A second 
technique couples the output of an OMT directly to two square-law detectors (Netterfield 
et al. 1995; Wollack et al. 1997). The beam is switched on the sky to measure the spatial 
anisotropy in two orthogonal polarizations. This approach measures the anisotropy in the 
Q Stokes parameters of the incident radiation field, and provides the most stringent upper 
limits on the spatial anisotropy of the polarization of the CMB. 

An alternate approach, employed in POLAR, is the correlation radiometer (Fujimoto 
1964; Rohlfs 1990). In this instrument the two polarization components are amplified 
in separate parallel amplifier chains; the output signals are correlated, resulting in an 
IF signal proportional to the U Stokes parameter. This type of instrument effectively 
"chops" between the two input RF signals at a frequency which is comparable to that of 
the RF signals themselves. The correlation polarimeter has a \/2 noise advantage over the 
Dicke-switched approach. An advantage of this differencing mechanism is that it has no 
magnetic or moving parts which have traditionally complicated experiments of this type. 
After a 45° rotation the correlator gives the Q parameter. POLAR rotates about the 
vertical in 45° steps at a few rpm. The rotation modulates the output sinusoidally between 
U and Q at twice the rotation frequency and allows the removal of an instrumental offset 
and other instrumental effects that are not modulated at this frequency. 

The sensitivity of the polarimeter is determined by low-noise HEMT amplifiers cooled 
to 15K by a commercial cryocooler. State-of-the art versions of these amplifiers achieve 
noise temperatures of ~ 10K over a bandwidth of ~ 10 GHz in both the K a and Q 
bands (Pospieszalski 1992; Pospieszalski 1995). This noise temperature is comparable to 
the antenna temperature of the atmosphere at a good observing site. Nevertheless, long 
integration periods are required to reach a sensitivity level ~ lfiK. The rms noise in a 
measurement of either Q or U (in antenna temperature) is given by the radiometer equation 
(Krauss, 1982). For the Q Stokes parameter: 



where T rec and T atm are the receiver and atmospheric noise temperatures, respectively, 
r is the total time spent observing the CMB; the time spent either Q or U is r/2. Av 
is the RF bandwidth and k = a/2 for a correlation radiometer. For the K a polarimeter, 
T rec ~ 20K and T atm — 10K, and Av is 10 GHz, resulting in a sensitivity to Q or 



K(T rec + T a t m + 2.7) 



(38) 
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U of NET = AQ rms ^Jr/2 ~ 460/zKs 1 / 2 . For the total polarized intensity we have: 
Ipo\ = VQ 2 + U' 2 . The error in 7 po i is A/ po i = \/2AQ rms = 650/iKs 1 / 2 , before foreground 
subtraction. To reach a signal level of A/ po i = lfiK for a single pixel requires an integration 
time of ~ 120 hours. Consequently, to measure polarization at a signal-to-noise ratio of 1 
per pixel for all 36 pixels, demands a total observation time of ~ 180 days. The Q-band 
radiometer requires a similar amount of time. 

3.2. Systematic Effects 

Because the anticipated polarization signal is a factor of ~ 10 times smaller than the 
temperature anisotropy currently being detected, the understanding of systematic errors is 
crucial. Polarization experiments have several advantages, however, that promise to make 
this effort possible. First, the atmosphere is known to be polarized only at a very low level, 
far below the expected level of CMB polarization (see section 4). Additionally, POLAR 
is essentially a total (polarized) power radiometer, which eliminates the comparison of 
pixels through different airmasses, and at different times. In anisotropy observations, 
beam switching often adds noise and additional chop-dependent signals. Potentially, 
atmospheric effects will have a smaller contribution to this type of experiment than to 
ground- based CMB anisotropy experiments and will allow longer observation times than 
have been possible in the past. Long-term observations are key to understanding and 
removing systematic effects (Wilkinson 1995; Kogut et al. 1996b; Bennett et al. 1993). 
Many spurious instrumental effects can be isolated from astrophysical effects by long-term 
integration tests with the horn antenna replaced by a cold termination. 

In Table 2 we list some important systematic effects encountered in previous 
polarization measurements and summarize the solution adopted by POLAR. A full analysis 
of all potential effects is obviously beyond the scope of this paper. In section 4 we discuss in 
detail the discrimination against foreground sources, such as: extragalactic sources, galactic 
sources, and the atmosphere. 

4. Foreground Sources 

A fundamental question for any measurement of the polarization of the CMB is 
whether the expected signal can be distinguished from polarized foreground sources. While 
astrophysical (non-cosmological) sources of polarized radiation are of interest for other 
fields, the measurement of CMB polarization is our main objective, so these sources are 
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spurious effects. These foreground sources all have spectra that are distinct from that of the 
CMB, and in principle can be distinguished from it by multi-frequency measurements. This 
technique has been employed for observations of CMB anisotropy (Brandt et al. 1994). 
However, polarized foreground spectra have not been studied as extensively. To estimate 
the intensity and spectra of these foreground sources, we rely on theoretical predictions 
and extrapolations from measurements at different frequencies of the antenna temperatures 
of these foregrounds. Here we summarize the properties of atmospheric and astrophysical 
(though non-cosmological) foreground sources. 

Synchrotron Emission. Diffuse galactic synchrotron radiation arises from ionized 
regions of our galaxy that posses magnetic fields. The antenna temperature of synchrotron 
emission obeys a power law: 

I synchrotron (^) ^ ^ j 

where (3 is referred to as the synchrotron spectral index. The polarization level n of 
synchrotron radiation is related to the spectral index (Cortiglioni & Spoelstra 1995): 

n 3/3 + 3 
3/3+1 

. Faraday rotation and non-uniform magnetic fields will reduce the level of polarization 
given by this equation. Full-sky polarization maps made at radio frequencies (Brouw 
and Spoelstra 1976) have been extrapolated to millimeter-wave frequencies assuming a 
power-law spectrum (Lubin & Smoot 1981). The radiation is linearly polarized between 
approximately 10% and 75%, depending on galactic coordinates. Below 80 GHz the 
polarized synchrotron emission dominates all sources, including the CMB if it is polarized 
at the 1 x 10~ 6 level, as shown in Figure |j. We estimate the spectrum of the synchrotron 
radiation by extrapolating the Brouw & Spoelstra (1976) measurement at 1411 MHz to 
millimeter wavelengths with the modified power-law spectrum used to fit the COBE DMR 
data (Bennett et al. 1992). For our modeling purposes we choose: —2.9 > (3 > —3.2 total 
intensity of approximately 50 /iK at 30 GHz, typical of high galactic latitudes (Kogut et al. 
1996a; Bennett et al. 1992), and n= 75%. . 

Bremsstrahlung Emission. Bremsstrahlung, or free-free, emission from ionized hydrogen 
(HII) regions is not polarized (Rybicki and Lightman 1979). However, bremsstrahlung 
emission will be polarized via Thomson scattering by the electrons in the H II region 
itself. The rescattered radiation will be polarized tangentially to the edges of the cloud, 
at a maximum level of approximately 10% for an optically thick cloud. The locations and 
emissivities of galactic H II regions are not well known, but Bennett et al. (1992) model the 
bremsstrahlung emission in the galaxy by subtracting a synchrotron model from microwave 
sky maps. In any case the polarization in the rescattered bremsstrahlung emission will be 



at least an order of magnitude smaller than the polarized synchrotron signal at frequencies 
greater than 10 GHz. We quote the result of Bennett et al. (1992) that 

T ,/-2- 15 

J - bremsstrahlung u i 

with total intensity ~ 40/iK at 30 GHz (Kogut et al. 1996a; Bennett et al. 1992). 

Dust Emission. The polarization level of interstellar dust is not well known. At 
low galactic latitudes thermal emission from dust particles dominates the near infrared 
spectrum. Depending on the shape and the alignment of dust particles, emission from dust 
particles may be highly polarized (Wright 1987). Using the dust spectrum measured by the 
COBE FIRAS (Wright et al. 1991) normalized to the IRAS 100 micron map, we find that 
high galactic latitude dust emission is negligible below 80 GHz, even when assumed to be 
100% polarized. We use the two temperature dust model (Wright et al. 1991) : 

At high galactic latitudes T dust ~ 10//K at 200 GHz (Kogut et al. 1996a; Bennett et al. 
1992). 

Extragalactic Point Sources. The dominant radiation mechanism for extragalactic 
radio sources is synchrotron emission (Saikia & Salter 1988). These sources have a net 
polarization of < 20%. Calculations made by Franceschini et al. (1989) of the temperature 
fluctuations in measurements of anisotropy of the CMB arising from unresolved, randomly 
distributed sources show that they contribute negligibly at 30 GHz to a 7° anisotropy 
experiment. If the orientations of the polarization vectors of these sources are uncorrelated 
over 7° regions, we would also expect a negligible contribution to the signal observed 
POLAR. We ignore the contribution of these sources in our foreground modeling. 

Atmospheric emission. The antenna temperature of the earth's atmosphere between 10 
and 60 GHz is dominated by an emission feature at ~ 22 GHz caused by atmospheric water 
vapor, and a series of emission lines at ~ 60 GHz due to molecular oxygen. In the absence 
of external fields, neither of these atmospheric components is known to emit polarized 
radiation in the frequency range of interest. However, Zeeman splitting of the energy levels 
of atmospheric molecules by the magnetic field of the earth can produce polarized emission. 
The valence band of water is completely full, and thus, does not exhibit Zeeman splitting. 
However, the 2 molecule has a non-zero magnetic moment due to its two unpaired valence 
electrons which interact with the Earth's magnetic field. We here discuss polarized emission 
from mesospheric oxygen, and show that it is negligible in comparison with the expected 
polarized intensity of the CMB. 
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The Zeeman effect breaks the energy degeneracy of the two unpaired valence electrons 
of molecular oxygen. The total angular momentum quantum number of the oxygen molecule 
is j = 1, which implies that the oxygen molecule's rotational spectral lines are Zeeman 
split into 2j + 1 = 3 distinct lines. Dipole radiation selection rules for transitions between 
these levels permit transitions as long as the change in magnetic quantum number, m, is: 
Am = 0, ±1. Transitions with Am = +1, for example, correspond to the absorption of 
a right circularly polarized photon or the emission of a left circularly polarized photon. 
The absorption and emission properties depend, therefore, on both the frequency and 
polarization of the radiation. The frequency of each Zeeman split level is (Liebe 1981): 



u z = v + 2.803 x 10" 3 J Br/(Am)[GHz] (41) 

where uq is the unperturbed frequency, rj is a shift factor with \rj\ < 1, and B is the 
magnitude of the earth's magnetic field, typically 0.5 Gauss throughout the mesosphere. 
The largest possible frequency shifts occur for r\ — ±1, which imply that the center 
frequencies for the polarized emission components will be confined to within 1.4 MHz of 
the unsplit center frequency. In principle, emission at these split frequencies could be 
up to 100% circularly polarized. Away from the center frequencies, the total intensity of 
emitted radiation decays with frequency as: / ~ (u-u ) 2 (R° sema, anz 1994). For a small 
shift in frequency, Au , away from the center frequency, the first order fractional change in 
emissivity can be shown to be: 

AI 2Az/ 
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For a single Zeeman split component, =1 = v _ z ^ m , where Av z ^ m = v z — Vq = 
2.803 x 10~ 3 -B?7(Am)[GHz], from equation [IT} To obtain the total contribution to the 



emission of both polarization components we must sum over left-handed and right-handed 
contributions: 

Aitot __ \p 2Av z ,Am 

I ~ a ~, , V-V Z 
Am=±l ^ 

However, we have for the shift factor in equation ^TJ: i](Am = +1) = —i](Am = —1) 
so the net effect on the emissivity is exactly canceled out by the two circularly polarized 
components (Rosenkranz 1994). Any second order contributions to the emission scale as 
~ (jj2 Vz f. i which implies a contribution of < 10~ 8 K for 26 < v < 46 GHz, i.e. the frequency 
band which POLAR will probe (Crill 1995). For these frequencies of observation there is 
also a small Faraday rotation of the plane of polarization of the CMB. Rosenkranz & Staelin 
(1988) show that the rotation of the plane of polarization will be less than ~ 10~ 2 degrees 
for these frequencies. Therefore, both the polarized emission and Faraday rotation of the 
atmosphere are negligible effects in the range of frequencies which POLAR will probe. 
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Of all the relevant foreground sources, only diffuse galactic synchrotron radiation and 
dust are expected to appear at a level comparable to the anticipated polarized CMB signals, 
(see Figure [|). In the next section we discuss techniques to remove spurious foreground 
sources using multi-frequency observations. 



5. Foreground Removal 

Here, we estimate our ability to subtract foreground sources in the presence of 
atmospheric and instrumental noise. We apply two different approaches to foreground 
removal. In the first approach we compute the anticipated error in our recovery of Q and 
U using an analytic technique developed by Dodelson (1995). In the second, Monte Carlo, 
approach we create simulated data and use a least-square fitting procedure to recover the 
underlying CMB Q and U parameters. We evaluate the quality of the fit by comparing the 
recovered and true CMB values. 



5.1. Foreground Modeling 

If there are several sources of polarized radiation at a given frequency, in a given pixel, 
the Stokes parameters describing the total radiation is the sum of the Stokes parameters of 
each source. Following Dodelson (1995), we represent the total signal as: 

3 3 

Q = Y, ? ipl + N 9 and U = 5] uT 1 + N n (44) 

i=0 i=0 

where the labels 0, 1, 2, 3 represent the CMB, synchrotron radiation, bremsstrahlung, and 
dust emission contributions to the total signal, respectively. The Q and U are expressed 
in terms of antenna temperature. The q l and u l are the thermodynamic temperatures of 
each signal source, and the F* are the spectral shapes of each source, following section 
[|. Q and U are written as vectors with dimensions equal to the number of observation 
frequencies. The components of Q and U are the Stokes parameters at that frequency, e.g.: 
Q = [Q(yi), QW, ■ " ■ , Q{vN obs )\- At frequency v. 
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(45) 

A similar equation holds for U. B V (T) is a blackbody spectrum with thermodynamic 
temperature T. N 9 and N u each represent the combined contributions from instrument 
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and atmospheric noise. For simplicity, we assume that the noise has the same Gaussian 
distribution for both Q and U, and is the same for each frequency channel. Furthermore, 
(N) = and (N 2 ) = a 2 . 

The unknowns in the signal are the coefficients q l and u l and the synchrotron spectral 
index, f3. To simplify our notation, we write these intensity coefficients as vectors: 
q = [q°, q 1 , q 2 , q 3 ]. Note that q° is the thermodynamic temperature corresponding to the Q 
Stokes parameter of the CMB, and is frequency independent. 



5.2. Error Analysis: Analytic Method 

Our goal is to recover q and u, particularly q° and u° from measured values of the 
Q and U Stokes parameters for each pixel at multiple observation frequencies. Given Q 
and U it is convenient to consider 7 po i, the linearly polarized intensity of the radiation 
Jpoi = \JQ 2 + U 2 , and the polar angle \a\ = |tan _1 Following Dodelson (1995), we 
wish to minimize the difference between the true CMB polarization signal and our best 
estimate of it. The fitting process is a x 2 minimization for the linear unknowns which 
produces our best estimates for the intensity coefficients: 

q' = n~^Tjq and u' = K~ lr q u 

where 

F*-F J ' , i Q F* i V-F i 



k ij = — and rf = — , rf u 
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The F* encode only the spectral shape (frequency dependence) of each foreground 
component. They are defined to be unit vectors such that F* • F J = 1 implies that F l and 
F J have the same frequency dependence. Note that the dimension of q' and u' may be 
different from that of q and u, if we choose to restrict our attention to fewer foregrounds. 
The upper limit on the dimension of q' and u' is the number of frequencies of observation. 

We define 

*U = VV) 2 + (u H ) 2 and \a H \= ^ 

as the ith component's polarized intensity coefficient and polar angle. The error in our 
analytic fitting process is: 

error(^ ol ) = y^" 1 )" = ^ (46) 

and 

error(|a'i) = - ^ (47) 
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We are interested in the recovery of the true CMB polarization intensity and orientation 
angle in the presence of the foregrounds. The error in the fitted CMB polarization intensity 
coefficient, error(ip 0l ), is found from equation pip to be the standard deviation of the 
system (instrumental + atmospheric) noise multiplied by a factor, e°, which depends only 
on the frequencies of observation and the choice of foregrounds for which one fits. It can be 
normalized such that its minimum value (i.e., with no foregrounds or system noise) is 1.0. 
The normalized parameter, known as the Foreground Degradation Factor (Dodelson 1995), 
is 

FDF = eVF° • F°. 

An optimized experiment will minimize not only the system noise a, but also the FDF. 
Another contribution to the analytic error defined by Dodelson (1995) is the error arising 
from the uncertainty in the spectral shape functions, F\ of the foregrounds, which vary 
with position on the sky. The situation is now more complicated because we must compute 
these contributions separately for both Q and U. This error, called o" s h ap e by Dodelson 
(1995), is negligible for our most important foreground, synchrotron radiation, even though 
we only know ft to 10%. 

5.3. Error Analysis: Monte Carlo Method 

To corroborate the results of the analytic error calculation, we perform a more explicit 
foreground removal simulation. The approach is similar to that of Brandt et al. (1994). 
First we choose a particular set of observing frequencies. Using equation fl4"5f) , at each 
pixel we create simulated signals for each frequency. Using foregrounds levels typical of 
high-galactic latitude regions. While the foreground temperatures are fairly well known, 
the polarization levels are not and we choose the most conservative estimates: 75% for 
synchrotron, 10% for bremsstrahlung, and 100% for dust. The CMB signal is chosen 
nominally at 10% of the anisotropy level (~ 3/iK), but the results of the simulation are 
independent of the exact value. For each of these signals we keep the true intensity, = i po i l , 
fixed but allow the orientation, = a 1 , to vary randomly. Noise is chosen from a gaussian 
distribution with standard deviation a in each frequency channel and added as in equation 
(fl5|). We then perform a least-squares fit to a particular signal and foreground model and 
determine our best-fit values for the recovered CMB intensity i po i /0 and orientation angle 
a' . We repeat this process 200 times, generating new values for the system noise and the 
a 1 each time. Finally, we compute the rms of the difference between the 200 input CMB 
coefficients, i po i°, and the corresponding best-fit CMB values, i® ol . A similar procedure is 
conducted for the true CMB polarization orientation angle, a , and the recovered angle, 
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a'°. In this process we average over many different relative orientations of the polarized 
foreground vectors. 

Typically the number of foregrounds (e.g., 4) simulated is greater than the number of 
frequency channels (e.g., 2 or 3), so we can fit for only some of the signals (e.g., the CMB 
and one or two foregrounds). As shown below, for a judicious choice of the observation 
frequency channels, the neglected foregrounds will not contribute significantly to the error 
in the recovery of the CMB signal. 

The primary unknown which is not determined by this analysis is the synchrotron 
spectral index, which varies with sky position. We estimate an uncertainty in the 
synchrotron spectral index by using a range of values for the fitting function F 1 ^) = v 13 . 
We evaluate of F 1 ^) for three values of the spectral index: (3 = —2.9, —3.05, —3.2. 

5.4. Results of Error Analysis 

Here we concentrate on a system noise range of 0.2 to 50 /iK (after a long integration 
on each pixel), and show the results of a two channel configuration observing at 30 and 40 
GHz, which is appropriate for POLAR. Here, the only foreground considered is synchrotron 
radiation, and we find that the degradation to our experimental sensitivity is: FDF= 2.7. 
Other frequency channel configurations show, very generally, that in order to fit for an 
extra foreground source (e.g. dust emission) without a severe increase in the FDF, one 
must observe more than three frequency channels. 

The plots of polarized intensity coefficient error versus system noise demonstrate that 
the quality of the recovery is fairly insensitive to the value of the input synchrotron spectral 
index. The Monte Carlo results are plotted as points, and the analytic result is plotted as 
a line. The excellent agreement of our analytic and Monte Carlo results confirms the idea 
that a successful observation should have a small e°, or equivalently, a low FDF. 

For CMB polar angle recovery, the plotted symbols represent Monte Carlo simulations 
at the two extreme values of synchrotron spectral index. Generally, the polar angle is more 
difficult to determine than the total linear polarized intensity, even at very low values of 
system noise. However, frequency configurations which yield an accurate recovery of the 
total polarized intensity will often give a more accurate recovery of the polar angle as well. 

From our analysis and simulation of a CMB polarization, we conclude that for a dual 
frequency experiment observing frequencies below 60 GHz, the most important challenge 
is to reduce the system noise. The contribution of the system noise is far more important 
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than that of combining bremsstrahlung and synchrotron in the fitting process or reducing 
the uncertainty in the synchrotron spectral index used to fit the data. We conclude that 
an experiment of the 30-40 GHz type with system noise lower than 1 fiK is capable of 
discriminating a ~ 1 — 3fiK CMB polarization signal from polarized galactic synchrotron 
radiation. 



6. Observation Strategy 

Constraints on, or detection of, the polarization of the CMB and its associated 
power spectrum depend greatly on the amount of sky coverage of the observation and 
the sensitivity of the radiometer. Sensitivity considerations are common to all CMB 
observations: time limitations restrict signal integration, and constrain the amount of sky 
coverage. We must reach a compromise between the integration time required to achieve 
the desired signal-to-noise ratio while also sampling a representative distribution of celestial 
regions. We now discuss our observing strategy in the context of the achievable level of 
sensitivity of POLAR. 

6.1. Sky Coverage 

The trade-off between the signal-to-noise ratio per pixel, s:n, and the total number 
of pixels arises frequently in designing CMB observing strategies. The ultimate goal is 
to discriminate between two hypotheses: (H and H po \), which are the null hypothesis of 
an unpolarized CMB, and the hypothesis that the CMB is polarized at a particular level. 
This discrimination is quantified by the confidence level and power of the measurement. 
POLAR seeks to make a primary detection, namely, a detection of a signal for which no 
prior detections have been made. There are four possible outcomes: the first two are 
correct detections of either Hq or H^. The other two outcomes are erroneous detections: 
a false positive or a false negative detection. Maximizing the confidence level minimizes 
the probability of a false detection, while maximizing the power minimizes the probability 
of making a false negative detection. Power quantifies the ability of an experiment 
to distinguish between competing hypotheses, and confidence quantifies the certainty 
of the detection. Traditionally, most CMB experiments have quoted only confidence 
intervals. Keating & Polnarev (1997) argue that the confidence level should be equal to the 
discriminating power for a primary detection. 

For a fixed total amount of observing time we wish to find the number of pixels which 
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maximizes the confidence level and the power of the experiment. We estimate the polarized 
CMB intensity from models described above. While our s:n per pixel is less than one (see 
section |3|), we expect that our signal-to- noise ratio defined over the total observation time, 
denoted by R, will be R ~ 1. For R — 1 it is shown that the optimum number of pixels 
is ~ 10, with a very weak dependence of the optimum on N. In order to minimize the 
susceptibility to certain systematic effects, as discussed above, POLAR observes a constant 
declination which corresponds to the zenith. The number of independent pixels which will 
be observed is N = 360° cos 6*i at /©a, where #i at is the latitude of Madison, WI and 0a is 
the antenna full width at half maximum, FWHM. We find N ~ 36, though we are free 
to pixelize our data into bins which are smaller than this, for instance pixelizing at the 
'Gaussian width' of our antenna, <tb = FWH M/2y/li\.2. Observing fewer pixels would 
increase the per-pixel s:n, but would either require tracking individual pixels over large 
angles on the sky or tilting the radiometer toward the north celestial pole. Either approach 
would undoubtedly introduce gravitationally modulated systematic effects into our data. 
Additionally, with TV ~ 36 we will have good coverage of the galaxy, which should allow 
for both the removal of foregrounds and a comparison with previous galactic polarization 
surveys (e.g., Brouw & Spoelstra 1976). For iV = 36 and R — 1, POLAR can expect to 
make a detection at the ~ 55% level of confidence and power for a CMB which is polarized 
at ~ lyuK. Of course, it is still possible to quote results with confidence arbitrarily close 
to 100%, but this would be at the expense of ability to reject the null- hypothesis (i.e., 
low power). In the context of a primary detection, requiring equal levels of confidence and 
power seems the most reasonable compromise. 

6.2. Sensitivity to the Power Spectrum 

The observation's sensitivity to CMB fluctuations is quantified by its window function, 
Wg. The observed two-point correlation function is related to the power spectrum and 
window function as follows: 

1 oo 

(Q(ni)Q(n 2 ) + C/(ni)C/(n 2 )) = ^-J2( 2£ + 1 ) x C f W ! x ^(cos0), (48) 

47r 1=0 

where, for example, Q(n) is the Stokes parameter measured for a pixel located in the 
direction n. Cj 1 is the power spectrum describing the degree of polarization on angular scales 
characterized by £, Wf is the window function of this observing scheme, cos(#) = ni ■ fi2 is 
the separation between pixels in this observing scheme. 

The analysis differs from that of most anisotropy experiments in several respects. 
The primary difference is that the observations are total-power in nature, rather than 
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differential. The window functions for this experiment will reflect the fact that there is no 
"chopping" of the beam in sky-position inherent in the observation. Single pixels will be 
formed by binning the acquired data, and differencing between pixels can be performed 
during analysis of the data; not during acquisition. This approach avoids systematic effects 
which can arise from mechanical chopping mechanisms. Data from POLAR will be analyzed 
using a variety of synthesized window functions, each sensitive to a different angular scale. 
In this respect the analysis will be similar to that of the Saskatoon Big Plate observations 
(Netterfield et al. 1995, Wollack et al. 1997). 

Window functions for observations with less than full-sky coverage are specified by 
three functions: the beam profile function, the beam position function, and the weighting 
or 'lock-in' function (White & Srednicki 1995). The beam profile function, G{9, O^as), 
where a B is the beamwidth, quantifies the directional response of the antenna, which is 
assumed to be gaussian. G(9, 9 i: a B ) effectively samples all angular scales larger than, 
approximately, the angular size of the beam. The angular coordinates of the center of the 
beam are specified by the beam position function, 9{, and the lock-in function, wf, is the 
weighting of each of the N binned pixels indexed by i, for the scan strategy denoted by a. 



We have G(9, 9 i} a B ) = -V exp ) , where a B = FWEM/2^2W2 = 0.052. 



where Pi are the Legendre polynomials, and if a (x = J2i w iG(9,9i,<j B ) quantifies the 
response of the antenna (for a differencing strategy indexed by a), when pointed in the 

direction of x. For a gaussian beam: H a (x) = J2i wf-^r ex P ( ~ ) • 



By varying the weight function we will obtain window functions "tuned" to sample 
specific multipole- space regions. For example, wf = (— 1)* +1 , differences pairs of nearest 
neighbor pixels. Pairs of pixels will be separated by a constant angle for each differencing 
strategy indexed by a, which runs from 1 to 18, corresponding to the N/2 distinct two-pixel 
difference window functions of POLAR. The number of unique two-pixel combinations, k, 
is given by k = ^ N ^y. m \ ' wnere N is the number of binned pixels, and m is the number of 
beams which are differenced. For N = 36 pixels, differenced pairwise (m — 2), k — 630. We 
note that the total power nature of POLAR, as well as the insensitivity of our experiment to 
atmospheric emission and thermal gradients, will allow us to perform two-pixel differencing, 
in contrast to the majority of ground based CMB anisotropy observations, which typically 
difference three or more pixels. 




Following (White & Srednicki 1995), the window functions are: 




(49) 




To estimate the rms polarization for the FWHM = 7° polarimeter with 36 pixels, we 
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have calculated the single pixel window function of POLAR, W° a . This contains only the 
diagonal elements of the window function matrix from (f49"l). This quantity will allow us to 
determine the relationship between the measured pixels and the underlying spectrum which 
is responsible for a particular realization on the sky. The polarization two-point covariance 
matrix (at zero-lag) for a particular theoretical model is given by: 

Crf 4^+!)* (50) 

To compute the theoretical rms amplitude, we extract the square root of the diagonal 
elements of C^f and obtain: 



n 



11118 



\ 



i-^(2f+l)C?W,< 



where the Cj 1 have been introduced in section |2|. We show this window function in Figure 
H]. Of course, the off-diagonal components of the window function matrix will contain 
additional information about cross-correlation between pixels, as opposed to W® a , which 
is the auto-correlation. The off-diagonal window function matrix elements will increase 
the effective signal-to-noise ratio of the experiment as a whole, and can be 'tuned' to 
sample specific £-space regions, up to the cutoff £ of the antenna, similar to the analysis of 
(Netterfield et al. 1995). For anisotropy experiments it is conventional also to quote the 
band power, which is independent of the details of the experiment. This approach facilitates 
comparisons between experiments, and contains information equivalent to that of the rms 
amplitude. 



6.3. Estimated Signal Level and Uncertainty 

The primary goal of POLAR is to measure the polarization of the CMB. We have shown 
in section 2, that the level of polarization is extremely sensitive to the ionization history 
of the universe, both before and after recombination. We expect, then, that the observed 
polarization signal will depend critically on the optical depth, r, for photons back to the 
last scattering surface. A preliminary estimate of the effect of reionization can be obtained 
by computing the expected rms polarization and associated experimental uncertainty for 
models of a reionized universe. We will now compute the effect of reionization on the power 
spectrum, and demonstrate that the characteristic signature of an early reionization is, in 
principle, detectable by POLAR. 

Figure |2] displays Cj 1 , for the power spectrum computed using CMBFAST (Seljak 
& Zaldarriaga 1996a) for various totally reionized (ionized fraction x = 1) scenarios, 
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parameterized by the redshift of reionization z r \. In Figure [7| we plot the expected rms 
polarization vs. z ri , for < z ri < 105, along with the statistical la uncertainties we expect 
based on our NET and observation time. The underlying power spectrum is a generic 
CDM model with Q = 1,Qb — 0.05, h = 0.65, A = 0, and pure scalar perturbations. 
The inclusion of a tensor component should enhance the large angular scale polarization 
(Crittenden, Davis, & Steinhardt 1993; Crittenden, Coulson, & Turok 1995), so this figure 
underestimates the rms polarization predicted by some cosmological models. This figure 
suggests that POLAR could begin to detect polarization of the CMB at the la level if the 
universe became completely reionized at a redshift z T \ > 45. 

7. Related Topics 

Although no detections of the polarization have been made, we can glean information 
about CMB polarization from anisotropy detections. In principle, these detections can be 
utilized to refine a polarization observing strategy. Here we briefly discuss the possibility 
of designing a polarization observing strategy, utilizing information from a well sampled 
anisotropy map. 

Anisotropy and polarization are invariably spatially correlated with one-another, and 
additional cosmological information may be obtained by studying polarization-polarization 
and polarization- anisotropy correlation functions. Correlation provides the only link 
between previous detections (of anisotropy), and the proposed measurements discussed 
in this article. CMB polarization can be decomposed into two components: one which 
is spatially correlated with the temperature anisotropy, and another, larger component 
which is uncorrelated. Ng & Ng (1996) and Crittenden, Coulson, & Turok (1995), 
demonstrate that, given a high-resolution CMB temperature map, it would be possible to 
identify celestial regions which are statistically more likely to posses higher levels of the 
correlated polarization component. As shown in Coulson, Crittenden, & Turok (1994), the 
uncorrelated polarization component dominates the correlated component by a factor of at 
least three. 

For detector-noise limited polarization experiments, it can be advantageous to search 
for polarization- anisotropy (QT) correlation in addition to polarization-polarization (QQ) 
cross- correlation. If the noise in the temperature anisotropy map is negligible in comparison 
with the noise of the polarization measurement, a, the error in (QT) will be linear in a 
while the variance in the polarization cross-correlation function grows as a 2 . In this limit it 
becomes advantageous to search for correlation. 
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Power spectrum generation programs, such as CMBFAST, compute (QT), along with 
the anisotropy and polarization spectra. This allows one to predict the distribution of 
polarization which is correlated with the anisotropy given a well-sampled anisotropy map. 
Finally, Crittenden, Coulson, & Turok (1995) describe potentially observable distributions 
of correlated polarization vectors on the sky which result from the velocity field of the 
photon-electron plasma and the type of metric perturbation which generates the anisotropy. 
If the angular resolution and experimental sensitivity of measurements improve by several 
orders of magnitude over the current levels, these patterns could prove to be direct 
observables of the microphysical properties of the photon-baryon fluid at the moment of 
decoupling. 

We also briefly discuss an effect which is peculiar to CMB polarization measurements: 
Faraday depolarization of the CMB by primordial magnetic fields. Faraday depolarization 
causes the plane of CMB polarization to rotate differentially due to a residual primordial 
magnetic field, which may have existed during the epoch of recombination (Basko & 
Polnarev 1980). The effect is akin to optical dichroism, familiar from the polarization of 
visible light. The net depolarization is frequency dependent, as the rotation of the plane of 
polarization of individual frequency components scales as ~ 1/V 2 . This effect is expected 
to be non-negligible below 30 GHz, for reasonable values of the primordial magnetic field 
(Harrari, Hayward, & Zaldarriaga 1997; Kosowsky & Loeb 1996). Detection of this effect 
is unlikely until the polarization of the CMB has been detected over a large frequency 
bandwidth, and with high angular resolution. Neither of the above mentioned phenomena 
are immediately relevant for polarization observations, but are nonetheless quite intriguing. 

8. Conclusion 

We have demonstrated that the detection of the polarization of the CMB is difficult but 
technologically feasible. A detection would permit the discrimination between heretofore 
degenerate theoretical predictions. Polarization of the CMB has a unique signature in 
both real and Fourier space, as well as distinct spectral characteristics. A detection of 
polarization, in conjunction with the current detections of CMB anisotropy could be the 
best available probe of the ionization history of the pre-galactic medium. This epoch 
of cosmic evolution is of great interest, and supplemental information from polarization 
detection could greatly advance our knowledge of the formation of structure in the early 
universe. The current generation of anisotropy measurements are sufficiently refined that the 
fundamental parameters of classical cosmology are beginning to be determined. Detection 
of polarization of the CMB also promises numerous dividends throughout cosmology, and 
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one readily observes that the status of polarization observations today is reminiscent of the 
status of anisotropy measurements a decade ago. 

The authors are indebted to many people who have worked on the planning and 
development of this measurement. Khurram Farooqui and Grant Wilson designed an early 
version of the instrument and observing scheme. Brendan Crill tested the polarimeter and 
calculated the atmospheric emission. Victor Derderian designed the rotation platform of 
the instrument. Melvin Phua and Nathan Stebor helped to diagnose systematic effects. 
Dave Wilkinson made key suggestions regarding the design of the correlation radiometer. 
Jeff Peterson reviewed the manuscript and made helpful suggestions. Conversations with 
Robert Brandenberger, Josh Gundersen, Ka Lok Ng, Lucio Piccirillo, Uros Seljak, and 
Matias Zaldarriaga have refined the focus of the project. This work is supported by NSF 
grant AST-9318785, a Ford Motor Company University Research Grant, a NASA GSRP 
Fellowship for BK, and a UTRA grant from Brown University for JS. 
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Fig. I. — Schematic Dependence of CMB Polarization on Angular Scale and Recombination 
Scenario. In non-standard models which predict a prolonged recombination, the angular 
scale and level of polarization are larger than for standard models. 

Fig. 2. — Polarization Power Spectra and Window function for POLAR. Three spectra are 
displayed for a CDM (Q — l,fl B — 0.05, H = 65,) universe with three different redshifts of 
reionization (z Ti = 0, 50, 100). Also shown (solid line) is the single beam window function of 
POLAR. 

Fig. 3. — POLAR Schematic. The correlation polarimeter is a wideband, total-power device. 
An orthomode transducer separates the two orthogonal linear polarization components of 
the incident radiation. HEMT amplifiers are cooled to 15K inside of a dewar. The remaining 
components are at ambient temperature. 

Fig. 4. — Polarized Foreground Spectra at Millimeter Wavelengths. Spectra of expected 
polarized radiation sources at high galactic latitudes, are shown for a 7° beam. A 3 /xK 
polarized CMB signal is shown, corresponding to 10% of the 10~ 5 CMB anisotropy. At 
frequencies lower than 90 GHz the polarization signal is dominated by galactic synchrotron 
emission (up to 75% polarized, as shown). Galactic bremsstrahlung radiation is not polarized 
in direct emission, but can be up to 10% polarized (as shown) after Thomson scattering. 
Galactic dust is shown conservatively with 100% polarization. 

Fig. 5. — Recovery of CMB Polarization Intensity Coefficient vs. System Noise. This plot 
shows the results of the simulations described in section 5. The vertical axis is the error in 
the recovered CMB polarization intensity for a single pixel, in /zK. The plotted symbols are 
the Monte Carlo simulation results for two different values of the synchrotron spectral index; 
the solid line is the analytic result for the Foreground Degradation Factor. The analytic 
and Monte Carlo results agree extremely well. This figure demonstrates the feasibility of 
detecting a CMB signal polarized at the ~ 1 — 3fiK level using a two-frequency configuration. 
For such a detection the system noise, defined over the entire observing time, must be below 
1/iK. 

Fig. 6. — Recovery of Polarization Orientation Angle vs. System Noise. This plot shows 
the results of the simulation described in section 5. The vertical axis is the error in the 
recovered CMB polar angle for a single pixel, in degrees. The plotted symbols are the Monte 
Carlo simulation results for two different values of the synchrotron spectral index. Similar 
to the polarized intensity recovery figure 5, accuracy in the recovery of the CMB polar angle 
depends primarily on the system noise, which is defined over the entire observing time. 
Recovery of the CMB polar angle with errors less than 30° is only possible if the system 
noise is lower than 1/iK. 
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Fig. 7. — Expected rms polarization in reionized universes. The solid line is the magnitude of 
the rms CMB polarization as a function of redshift of reionization. All reionization scenarios 
assume total ionization (x = 1), except for z ri = 0, which represents no reionization. The 
gray band represents the la experimental uncertainty for POLAR, observing 36 pixels for a 
total time of 1 x 10 7 seconds, with NET= 460/iKs^ after simulated foreground subtraction. 
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Table 1. Experimental Limits on Linear Polarization (95% Confidence Level) 



Reference 


Frequency (GHz) 


Sky Coverage 


Limit Tp i/T cmb 


Penzias & Wilson 1965 


4.0 


scattered 


0.1 


Cadcrni et al. 1978 


100-600 


near galactic center 


0.001-0.01 


Nanos 1979 


9.3 


declination= +40° 


6 x 10~ 4 


Lubin & Smoot 1979 


33 


dclinations 38°, 53°, 63° 


3 x 10~ 4 


Lubin & Smoot 1981 


33 


11 declinations -37°to +63° 


6 x 10~ 5 


Partridge et al. 1988 


5 


43' x 43' region, declination 80° 


4 x 10~ 5 


Wollack et al. 1993 


26-36 


Degree scales, about NCP 


9 x 10~ 6 


Netterfield et al. 1995 


26-46 


Degree scales, about NCP 


6 x 10~ 6 
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Table 2. Expected Systematic Effects 



Effect 



Mechanical Strain a 
Magnetic Coupling 15 
Microphonics 

Electromagnetic Interference 13 
Offsets 

Radio-Frequency Interference* 
Thermal Variations 5 
Sidclobc Pickup 11 



Origin 

Instrument Rotation 
Rotation in Earth's Field 
Mechanical Vibration 
Local Sources 
Polarization Cross-Coupling 
In Band Sources 
Diurnal /Environment 
Sun/Moon/Earth 



Removal / Control Method 

Vertical Drift Scan 
No Ferritc Components 
Mechanical Isolation 
Shield/Filter 
Isolation/Instrument Rotation 
Identification and Data Editing 
Temp Control/ Shielding 
Low Sidelobe Antenna/Shielding 



a A problem with any radiometer that must move in a gravitational field is position-dependent stress and strain on 
waveguide joints, etc. In POLAR these problems are minimized by staring at the zenith, so no gravitational torques 
are present. The rotation speed is slow, ~ 1 rpm so accelerations on stopping and starting rotation are small. 

b A particular concern is the coupling of the Earth's magnetic field to the radiometer. The COBE DMR had 
ferritc Dickc switches which produced a spurious signal at the ~0.1 mK level (Kogut ct al. 1996b). POLAR has no 
ferrite components, but other components such as amplifiers, etc., may have a low-level magnetic field dependence. 
Modulation of these effects can be minimized by maintaining a constant orientation of rotation axis with respect to 
the Earth's field 

c The effects of vibrations that occur during rotation are reduced by taking data while the instrument is stationary 
and by stiffening the support structures. 

^This effect can be controlled by Faraday shielding the instrument and by filtering electrical lines. 

c Cross polarization in the antenna and the OMT create correlated signals in both arms of the polarimeter, which 
in turn produce an offset signal. If this offset is stable on the timescale of instrument rotation it is subtracted by 
phase-sensitive detection at the rotation frequency. 

c 

RF sources that occur in the radiometer RF band or IF band arc becoming increasingly troublesome. Of particular 
concern in the future will be communications satellites operating in the bands of interest. 

s Temperature variations in the radiometer that occur at the rotation rate of the instrument can be mitigated by 
active temperature control and by shielding the instrument from the Sun. The latter function is naturally performed 
by the ground shields so that the antenna and receiver are completely shielded. 

''The polarimeter must be able to reject or discriminate against emission from the Sun, Moon, and Earth, which 
appear only in the sidelobes of the beam. None of these sources are expected to be significantly polarized, but 
asymmetry in the antenna response to the two linear polarizations will create spurious signals. Requiring the total 
power from these sources to lie below 1 /iK demands 75 dB and 63 dB sidclobc rejection for the Sun and Moon, 
respectively. Assuming 20 dB rejection from the ground screen, this level of rejection can be achieved with the 
corrugated horn antenna (Janssen et al. 1979) if data are rejected when these sources lie closer to the zenith than 
50° and 30° respectively. Binning of the data in Sun-centered or Moon-centered coordinates will allow us to uncover 
correlations between the position of these objects and the response of the polarimeter. 



VACUUM 



WINDOW 




Microwave Polarization Signal Sources 




Frequency in GHz 



30 and 40 GHz 

5Q r~ i i i i i i i i | i i i i i i i i |~~ 



40 - A A A 

: A H □ D □ 



fee 

T3 



i 20 



A 



A 



□ □ 



□ 



A : Synchrotron S.I. (8 = -2.9 



□ 
A 



10 — □ : Synchrotron S.I. (8 = -3.2 

FDF = 2.7 

L I I I I I I I I I I I I I I I I I I I I I I I I I 

0.1 1.0 10.0 100.0 

a [/xK] 



2.0 




